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STEADY  STATE  OF  A NONLINEAR  EVOLUTIONARY  EQUATION 

by 


Joel  C.  W.  Rogers 
Applied  Physics  Laboratory  of  the 
Johns  Hopkins  University 
Johns  Hopkins  Rd. , Laurel,  Md.  20810 


A lecture  presented  at  Institut  de  Recherche  d' In forma tique  et  d’Automatique 
in  Rocquencourt , France,  on  May  9,  1978.  This  lecture  is  included  in  the 
publication  Seminaires  IRIA  Analyse  et  Contr6le  de  Systemes  - 1978. 


1.  introduction 

in  a recent  paper  (Rsf.  2),  Brenis,  Berger,  and  Roger 3 have  presented  an 


algorithm  for  solving  the  initial  value  problem 

-t-  Lf(u)  = 0,  x e -Tl  , 0 < t«i  t,  (1.1a) 

3 f(u)  - 0,  x t 3il  , 0 ^ t ^ T,  (1.1b) 

u(o)-.  u0€  n L°°CJl)  , x&fl  , (1.1=) 

in  a bounded  domain  Cl  with  smooth  bo'undary.  f satisfies  f(0)  m.  0 and 

0 £ f (u)  - f (v)  £ u - v for  u v.  (1.2) 


L:  D(L)  C L'Cil)  — ^ L’(il)  ls  a c-osed  linear  operator  whose  domain  D(l)  is 

dense  in  iJ C Cl)  and  whoso  as-jociatcd  somi-group 

- Lrt 

s(t)  r e,  . t > 0,  (1.3) 

1.3  contractive  in  L'CCl)  and  L (Cl)  . B is  a linear  operator,  such  that 
the  solution  of 
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L h m g,  xtfl, 

(1.4a) 

3 h - 0,  x&  3/1  , 

(1.4b) 

satisfies,  for  all  g g.  R(L), 

llh"L’ui)  4 “•  IKllu'(n.) 

for  sctne  <*.  > 0. 

I 

(1.4c) 

The  version  of  the  algorithm  which  is  used  as  a basis  for  deriving  the 
principal  theorem  of  this  paper  is  as  follows.  Given 

u — uc  , (1.5a) 

one  constructs 

n n 

u ~ Fo  uol  n:  1,  2,  3,  (1.5b) 

where 

F0  u - f(u)+  S f(u).  (1.5c) 

(We  use  the  symbol  **3"  for  S(^)  when  there  is  no  chance  of  confusion.)  un 
approximates  u(n'7/ ) in  the  following  sense: 

Hu"  - > ° 33  r"  "rT  - ^ (1.6) 

uniformly  for  t in  the  interval  [0,t]  (Ref.  2). 

We  are  interested  in  the  inhomogeneous  equation 

Lf(u)  -g(t),  Jt-H.  , t r 0,  (1.7a) 

3 f(u)  - 0,  x fe  , t > 0,  (1.7b) 

u(°)  - u0  , x fe  il  , (1,7c) 

particularly  in  the  st3udy  state  of  (1.7),  when  one  exists.  Thus,  some  en- 
largement of  t.ue  results  given  above  is  in  order.  The  generalizations  required 
involve  the  extension  of  equation  (1.6)  to  the  inhomogeneous  problem  and 
replacement  of  the  upper  time  limit  t by  00  . jn  some  cases  it  may  also  be 
desirable  to  drop  the  restriction  that  Ci  be  bounded. 

(ff!  .Oil  il  13  unbounded,  we  add  to  (1.1)  and  (1.7)  the  asymptotic 
condition 

A f M °-  (1-9) 

lx  V — » 

where  a is  a linear  operator,  and  we  ass'un9  that  for  any  bounded  domains 


* 
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<=■  iX  and  Sc  n and  function  g € R(l)  with  supp  g c S'  , 
3 oc(  <8 ' ) > 0 such  that  the  solution  of  (1.4a,b)  and 


xcil 

(X  1 -*oo 


satisfies 


1 ^ L'(8)  ~ ( & > & ) H 8 H 


(1.9a) 


(1.9b)) 


Hera  we  will  not  give  any  formal  proof  that  the  algorithm  of  Ref.  2 can 
be  used  in  these  extended  circumstances,  instead,  we  shall  assume; 


(Al)  For  g £ R(L)  and  all 


S(t)  - g + g» (t) 


(1.10a) 


such  that 


jlu.ll 

JQ 


d't-  < o° 


J0MS  "L'Cil) 

(1.7)  has  a solution  which  approaches  the  steady  state  value  u satisfying 


(1.10b] 


L f (u)-n  g,  xeil  ; 

(A2)  uniformly  for  'l'  & (0,7^]  » the  quantity  un  computed  by 


u - u0  , 


rvn 

a - 


n r 

u 4-  j S(  Y - £ ) g»(  £ + nT)d£ 


where 


u - u - f(u)  + S( r ) * (4)  + g T 


(1.11) 


(1.12a) 

(1,12b) 


(1,12c) 


| f __ 

8*(T)  - y S(t)dt  g,  (1.12d) 

approaches  in  L‘  C/L)  , as  T'' n 4 , the  steady  state  u*('7/  ) 


satisfying 


l-S  cr) 


f(u*(r))  - 8* » 


(1.13) 


for  all  g which  satisfy  (1.10); 

(A3)  Given  any  T > 0,  (1.6)  holds  uniformly  for  t£  []o,t]  , where  un 
and  u refer  to  solutions  of  (1.12)  and  (1.7),  respectively. 

We  note  that  the  assumptions  (Al)  - (A3)  are  not  completely  independent. 
For  example,  (a2)  and  (a3)  obviously  imply  (Al).  in  addition,  one  may  derive 


from  (Al)  and  (A3)  a weakened  form  of  (a2)  which  is  expressed  as  follows; 


- 4 - 


(AS*)  For  any  £ > 0 ^ T > 0 and  T & (0,Yl  such  that,  for  the 


quantities  u computed  by  (1,12)  and  the  steady  state  u of  u given  by  (1.7), 


l!“"  -ullL'OL)  ^ 


(1.14) 


when  n Y ^ to  and  Y . 


We  will  see  in  the  sequel  that  (Al),  (A2»),  and  (a3)  (and  thus  (a1)  and 
(A3))  will  be  sufficient  for  us  to  prove  our  principal  results  pertaining  to 
the  steady  state  solution  of  (1.7).  to  derive  (A2*)  from  (Al)  and  (a3),  note 


first  that  (1.10b)  and  (1.11)  imply  there  is  a tq  such  that  for  T ^ Tc 

L lls,lli_!cil) 


(1.15a) 


iu(T)  - u ll  l'c  A.)  ~ 3 


(1.15b) 


next,  note  that  the  operator  F^  in  (1.12c)  is  LC-fl)  -oontractive : 

HF3  U " F3  vllL-ai)  - " vl(L'CA)  • f1*16) 

(1.16)  and  the  contract iveness  of  g(t)  in  (JOT)  , when  combined  with  the 

bound  (1.15a),  imply  that  the  functions  computed  by  (1.12b)  satisfy 

IU"  <*)  - u"*  11  ^ § (1.1V) 

for  n ^-no  and  n0Y  to  • From  (a3)  with  t replaced  by  T0  +■  , it 


follows  that  3 Tj  such  that  for  0 'Y  ^ TJ 

llun°  - u(noy  )llL.cil)  ^ | ( 1 . 18a ) 

where 

T0  £ no  r < T0  + ro  . (1.18b) 

Through  (1.18)  and  (1.15b)  we  bound  l(un°  - u 1 1 C ) * All  we  need  do 

is  extend  this  to  a bound  on  l|  F^«n°  uf1°  “ u^L!(.-fl-)  and  USe  t0 

get  the  desired  result  (1.14).  To  provide  the  missing  link,  we  observe  from 
(1.11)  and  (1.12d)  that 

I«*7':(I-3(T)!I=i(I-  3 (r)K(n), 

and  thus 

r-^f(u)ng*,  (1.19) 

on  account  of  (1.4c)  or  (1.9b).  Hence,  from  (1,12c), 


f9,  u-  u 


(1.20) 
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V"'  U"‘  -ullL'Ul) 


4 u"°  -u  „ 


ucn.) 


(1.21) 


follows  from  (1.16). 


/although  (a3)  Is  an  assumption,  for  Q-  bounded  with  smooth  boundary  and 
relatively  mild  conditions  on  g*,  it  is  only  a trivial  extension  of  the  conver- 
gence theorem  cited  above  for  the  homogeneous  problon  (1.1)  (Ref.  2).  yor 
example,  sufficient,  but  certainly  not  necessary,  conditions  on  gt  to  establish 
(A3)  would  be 

( 3 up  g»(t  +■  y)  - inf  g’(t+  y))dxdlt-»0  as  n— >0.(1, 22) 

-lo  Aft  °-  n / 

For  11  unbounded,  one  may  prove  (a3)  by  approximating  with  arbitrary  accuracy 
in  L1  C-H.)  * for  QQy  T > 0,  the  solutions  of  (1.12)  with  nft  T by  functions 

which  are  similarly  defined  in  terns  of  a boundary  value  problem  with  homoge- 

/VI 

neous  boundary  data  in  a suitably  large  bounded  subdomain  _Q  (t)  C.  _fL 

With  regard  to  the  assumptions  (Al)  or  (A2)  relating  to  the  approach  of 
solutions  of  (1.7)  or  (1.12)  to  a steady  state,  we  note  that  these  assumptions 
may  be  verified  directly  in  particular  cases  of  interest,  such  as  the  one-phase 
Stefan  problem  in  the  absence  of  sources.  More  generally,  when  the  operator 
3 (t ) is  strictly  contractive  in  L'(.H)  and  for  sane  a > 0 

l|s(t)llcai)  * (1*23) 

an  adaptation  of  the  argument  which  led  to  (1.16),  in  which  we  replace  y^  by 
Fg<*  , u by  un  , and  v by  u°  ' , enables  us  to  establish  (a2).  when  the 

assumption  of  strict  contractiveness  la  dropped,  one  may  still  derive  a partial 
result  regarding  the  approach  to  steady  state  of  the  quantities  generated  by 
the  algorithm  (1.12)  (cf.  Lemma  1 of  the  next  section). 

We  will  show  when  (Al)  - (A3)  hold  that,  under  suitable  further  restric- 
tions on  g* , the  steady  state  value  "u  of  the  solution  of  (1.7)  coincides  with 


the  steady  state  value  u of  the  solution  tj  of  the  problem 

n-t  + l t(u)  - g,  xfe  A , t 7 o, 

B f(\ 7)  - o,  x fc  9A  , t ■>  0, 

7(0)  - uQ  +-  \ g»(t)  d±  , x £ Q . 


(1.24a ) 
(1.24b) 
(1.24c) 
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A proof  which  is  general  with  respect  to  operators  i and  functions  f with  the 
properties  (1.2)  and  (1.3),  but  with  sons  restrictions  on  g*  and  u0  , will  be 
given  in  the  next  section,  our  method  will  be  to  establish  that  a similar 
result  holds  exactly  for  the  functions  generated  by  the  approximate  algorithm 
(1.12),  under  assumption  (a2),  and  that  such  a result  holds  approximately  under 
assumption  (A2’)«  In  either  case,  the  property  is  established  for  the  exact 
solution  by  using  the  convergence  of  the  approximate  solution  to  the  exact  one, 
as  assumed  in  (a3).  in  the  third  section  we  will  discuss  a direct  proof  of  the 
theorem,  without  reference  to  our  approximation  of  the  exact  solution  and  the 
convergence  of  the  approximation,  for  the  special  case  of  a one-phase  Stefan 
problem  with  g ~ o.  This  will  serve  to  illuminate  the  sorts  of  problems  we  may 
expect  in  trying  to  obtain  a direct  proof  for  the  more  general  case.  The  fourth 
section  produces,  for  the  case  of  a one-phase  Stefan  problem  with  gm  0,  first 
a counter-example  to  one  possible  generalization  of  the  theorem  of  section  2, 
and  then  a proof  of  a weakened  form  of  the  generalization.  The  final  section 
gives  an  application  of  the  results  obtained  for  the  one-phase  Stefan  problem 
with  g m 0,  whereby  it  is  shown  that  a time-dependent  free  boundary  problem 
which  arises  in  the  theory  of  anodic  smoothing  is  equivalent  to  the  steady 
state  of  a one-phase  Stefan  problem. 

Before  proceeding,  let  us  look  at  the  problems  (1.7)  and  (1.24)  whose 
steady  state  solutions  we  compare,  in  either  case,  we  get  equation  (1.11)  for 
the  steady  state.  Because  of  (1.4)  or  (1.9),  f(u)  is  thus  uniquely  determined. 
Accordingly,  if  f is  1-1,  all  our  results  reduce  to  a triviality.  The  only 
interesting  case  is  when  f(u),  constrained  to  be  monotone  by  (1.2),  is  constant 
over  an  interval  of  values  [u"  ,u+]  of  u.  stefa.n  problems  are  t*.e  best-known 

p rob lass  in  this  category.  Tith  h(x)  given  as  the  solution  of  (1.4a,b)  or 
(1.4n,b)  and  (1.9a),  we  will  find  it  convenient  to  define  the  functions  u-  (x) 
by 

u+  (x)  r.  sup[*  | f ( £ ) n h(x)^  , (1.25a) 

a-  (x)  inf  j ? | f(*  ) = *(*)  l • 

. n -» 


(1.25b) 


~ -nrJMSSTL:- 
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It  is  apparent  from  the  outline  given  above  that  at  the  current  time  our 
progress  regarding  the  steady  state  problem  has  been  uneven.  This  is  because 
the  results  reported  here  were  originally  obtained,  in  the  course  of  work:  on 
water  waves  (Ref.  9),  for  the  one-phase  Stefan  problem,  and  at  prerent  they 
have  been  extended  to  the  more  general  problem  (1.7)  only  in  part.  The  incom- 
plete state  of  this  work  is,  of  course,  also  reflected  in  the  fact  that  we  have 
found  it  necessary  to  assume  (Al)  - (A3)  in  order  to  obtain  the  principal  theo- 
rem of  this  paper,  relating  the  steady  states  u and  "5  of  (1,7)  and  (1.24* , in- 
stead of  obtaining  the  result  directly  without  such  assumptions. 

indeed,  the  restriction  of  uD  and  Jo  tgMdt  to  [_'  (X-)  is  probably 

unnecessary  when  XI  is  unbounded,  and  it  seems  quite  likely  to  us  that  assump- 
tion (a2)  can  be  dispensed  with  and  replaced  by  the  premise  of  the  following 
conjecture. 

Conjecture:  if,  for  all  u with  supp  u c S’  , & bounded  and  Sfc.  XI.  , 

/V  £ AJ 

there  is  a tfe  ( ,gr  , ft  ) such  that,  for  & bounded  and  S5  C XL  , 

ll5Wu||L'<m  ^ l|ullL'ce) 

whan  t te  ( £ , ff  ) , then 

Hmax(u0  - u+  ,0)  + max(u_  - ,°)((  ^ < c« 

implies  that 

Hinax{(pJ]  uo  ) - u+  ,0)  + mx(u"  - (y*J  uo)>0)lL  ~ -^*0  as  n— 

for  any  bounded  domain  £ c.  XI_  . 

At  present  we  have  not  satisfied  ourselves  as  to  the  validity  of  this  conjec- 
ture. 


As  I write  this  up,  it  occurs  to  me  that  the  question  of  the  steady  state 
of  (1,7)  really  involves  only  two  things  — the  linear  operator  i and  the  func- 
tion f — and  that  it  should  be  possible  to  derive  all  our  results  in  a more 
direct  and  simple  way,  in  a manner  which  can  be  extended  to  the  case  when  u is 
a vector,  solely  in  terms  of  properties  of  tne  linear  operator  and  the  non- 
linear algebraic  function.  However,  I leave  that  to  the  future. 


7 
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2.  The  gase  of  positive  sources" 

The  main  task  of  this  section  is  to  prove  the  following. 

Theorem  1:  under  assumptions  (Al)  - (A3)  or  (Al) , (A2»),  (A-3),  and  condi- 
tions (1.10)  and  (1.22)  on  g,  if  in  addition 

uQ  (x)  » u"  (x) , x e CL  , (2.1a) 

and 

g*(x,t)  ^ 0,  x fell  , t > 0,  (2.1b) 

the  3t9ady  state  u of  the  solution  of  (1.7)  is  the  same  as  the  steady  state  u’ 
of  the  solution  of  (1.24).  (u-  (x)  is  defined  by  (1.25b).) 

We  call  this  the  cuse  of  "positive  sources"  because  of  the  inequalities 
(2.1a)  and  (2.1b).  as  we  stated  in  the  introduction,  the  proof  of  the  theorem 
will  be  obtained  by  referring  to  the  algorithm  (1.121  and  examining  the  proper- 
ties of  its  steady  state,  or  approximate  steady  state,  according  to  whether  we 
assume  (A2)  or  (a2').  First,  we  give  a lemma  which  shows  that  a quantity  close- 
ly related  to  un  in  (1.12)  always  has  a steady  state. 


Lemma  1;  when  f satisfies  (1.2),  g*  is  given  by  (1.12d1  with  g£  F(l), 


and  F^  ia  given  by  (1.12c),  for  u £ L'C-TL)  the  qiantities 


In  - max(min(F.ni  u,u+  ) ,u“  ) 


approach  a limit 


'I  " V U*'  L'Cil) 


0 as  n — > 


(2.2a) 


(2.2b) 


Proof;  w®  use  the  L -contractiveness  of  F^»  , as  expressed  by  (1.16), 
and  the  fact  that  u £ f-'  (h)  where  h satisfies  (1.4e,b)  (and  (1.9a)  where 
applicable)  13  a fixed  point  of  f^»  . as  expressed  by  (1.201.  since  from 
(2.2a)  and  (1.85)  | £ f1  (h), 


V “ ■ *«  **  L’tn.)  - "V  “-inllLVU) 


f 2.3' 


ll;V  u-{-W  ll“I,r9‘ u,:i*  ’ "u\<iTlltt'  u,u_  ’"fca) 


lfV  u ""“'V  u'“r’  ‘ U-"AW 


I 
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Insarting  these  expressions  into  (2. 3)  and  repeating  the  result  for  n+-l,  n+2, 

• . • , n + m - 1,  we  obtain 
m 

1L  H^p  -ini-p-,  llmx{p  ,n+rT1  u,u+  ) - uMt^C-fT) 


• 

. . . , n + m - 1,  we  obtain 
m 

^ ^ “$n+p-| 

1(  max(F3' 

+ Hu"  -min(F9rm 

u,u- 

^ Hmx<rar  u,u^  ) - u+HL,tiI)+  \\ u-  - min(F3n,  u,u-  )||Ll(:il^  . • (2.4) 

(2.2b)  follows  froa  (2.4)  immediately. 

Remark  1;  (2.2b)  defines  the  operator  F * . 

Remark  2j  under  assumption  (A2),  the  quantities  F « u approach  t » u as 

3 3 

a °°  • Under  assumption  (a2»)  , for  6 > 0,  n r ^ To  ( €-  ) , and  'Y L (£), 

1 1 max(F^  u,u+  ) - ^+l^c_a)  + Hu"  - min(j£  u,u“  ) ||  ^ 6 - (2-5) 

Remark  3:  it  follows  from  (1.12c),  (1.2),  and  the  monotoniclty  of  s that 
F^»  is  a monotone  operator,  jf  u ^ u-  , the  ?n  given  by  (2.2a)  satisfy 

> (2.6) 

and  then  the  convergence  expressed  by  (2.2b)  holds  also  in  L^C-O-)  if 

u + € [f°(A) 

The  next  lenna  givs3  our  basic  result  for  positive  sources. 

Lemma  2:  when  f satisfies  (1.2),  g*  is  given  by  (1.2d)  with  g 4.  R(L) > F^ 
is  given  by  (1.12c),  and  f^*  is  defined  by  (2.2b),  we  huve , for  uo  > u~  and 
u,  ^ 0, 


11 V ,u“"  1 • V ‘V  u=  + ">  >>W) 

proof;  Let  us  define,  for  integers  n ^ 0, 


- 0. 


An'  V (U°  + U|  ): 

_ _ n 


An^  ~ V (J>  U°  + U) 

We  will  show  that,  for  p half  of  a non-negative  integer,  we  can  write 

(.p)  (.0) 

AP  ~ Uo  + u, 


(2.7) 

(2.3a) 

(2.3b) 

(2.9a) 


- V uo  + u, 


whore 


a0<|J)  ^ u“  , u((f>)  ^ 0. 

The  re3U.lt  is  obviously  true  for  p — 0 if  we  write 


(2.9b) 


(2.13) 
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1 


Co)  Co) 

u0  - Ue  , u,  m u,  . 


(2.U) 


If  we  can.  prove  that  there  is  a decomposition  (2.9),  (2.10)  for  p = '/z  , a 
similar  proof  will  establish  the  result  for  all  suitable  p. 


-i  c ) 

We  want  to  show  that  i u ^ u"  such  that 


Ai  “ A,/x  - ^ Uo 


A./z  > uo 


From  (1.12c), 


A,  - AVj_  - (S  - I)  f(uo  -t-  u(  ) - (3  - I)  f(u0  ) 


(2.12a) 


(2.12b) 


(2.13a) 


and  from  (1.12c),  (1.19),  (1.11),  (1.25),  and  (1.4a,b)  (and  (1.9a)  if  appropri- 


ate), 


C 'll.)  t'/i.)  (i/.) 

*3.  Uo  - Uo  - O - I)  (f(Uo  ) - f(u“  )). 


(2.13b) 


We  can  make  the  right-hand  sides  of  (2.12a)  and  (2.12b)  equal  if  we  choose 


“o/0  - lnf  [ rl  + U,  ) - f(uo  ) + f(u-  )]  . (2.14) 

Since  u,  ^ 0,  it  follows  from  (2.14)  and  (1.25b)  that  ^ u"  . to 

establish  (2.12b),  we  find 

f(A  t/J  - f(uo  + u(  + (S  - I)  (f(uQ  ) - f(u-  ))) 

> f(uc  + u,  - f(uo  ) *■  f(u"  )) 

^ f(uQ  + u,  ) - f(uQ  ) +■  f(u-)  n f(uoC',1-)),  (2.15) 

upon  use  of  the  contractiveness  of  s in  L^C-O.)  (monoto.nicity) , uo  ^ u"  , the 
property  (1.2)  of  f,  and  (2.14).  comparing  with  (2.14),  we  get  the  desired 
result.  Hence  we  may  regard  (2.9)  and  (2.10)  as  established. 


(2.15) 


From  Lemma  1,  in  L'C  fL)  min(An  ,u+  ) — =»  F0*(u0+-  u,  ) as  n— * ^ . 


Thus , 


(ninfF  , An  ,u+  ) - min(y  , min(An  ,u  + ),u+  jldx  — > 0 as  n— * . (2.16) 

JAn^+ 

Aa  was  pointed  out  in  Remark  3,  f„*  is  monotone,  so  is  the  operator 
min(Fg,  ( • ),u+  ).  Thus,  if  wn  is  a function  satisfying 

“in(An  ,u+  ) £ wn  4.  An  , (2.17a) 


it  follows  from  (2.15)  that 


* 

. wn<u+  w"  ) - wn]‘ 


0 as  n — » oc  , 


(2.17b) 
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or,  with  reference  to  (1.13c), 


Pick 


w. 


r 

[(3  - 1)  f (w„)  + 3"f] 

- = 

i 

Jwn<ur 

+ 

S"' 

C 

w O 

in) 

uv 

“.W’.  < 

u+ 

»o<m  < ut  ur 

U0l"> 

a3  n — > c*0 


(2.17c) 


u . 


i") 


^ u 


(2.18) 


It  is  clear  that  wn  satisfies  (2.17a).  in  addition 

tn) 


f(wn  ) = f(ue 


(2.19) 


on  account  of  (2.10).  Accordingly,  from  (2.17c), 


An<  u 


[(3-i)  f(u0w)  + g»r]ax: 


Where  An  = u0cr^  + u)cn)  "u+‘  , the  monotonicity  of  f^*  implies  that 
uo^  4-  u,cn)  > min(u0tr')  ,u+  ) +-  u,Crrt  5^  u*  , (2.21) 

or  Anvi/l_'3>  u+  i from  (2.9b).  integrating  this  result  over  the  set 
£ x | An  (x)  u"1"  ] and  adding  it  to  (2.20)  we  get,  on  referring  to  (2.9), 

I min.(A 0<-,/q_  >u+  ) - min(An  )|<ix  > 0asn->oo  . (2.22) 

Recalling  the  definitions  of  An  and  An+y  in  (2.8),  we  see  that  the  lenma  is 
proved. 

Upon  repeated  application  of  (2.7)  we  get,  when  uc  u~  and  u^  0, 

1 ^ i £ m, 


cn) 


iv-"1 

A„<u* 


- U0(n)  cix->0  as  n- 


oo 


(2.:- 


on  \ 


f 


v 


. m 

m-| 

UC+  Fg*  U 

+• 

1 

Um_(  4-  um) 

** 

3 

1 

1 u 4-  F 

0 x a* 

rn  3 

u« 

4-  ...  + F,  4-  Um_,  4-  Um) 

V ( 

Z u-  ) in 

L-O 

L1 

ca.) 

(2.23' 

To  prove  Theorem  1,  all  we  need  to  do  is  make  several  observations . 

. I 

First,  from  the  convergence  assumption  (a3)  and  the  L -contract lveness  of 

, it  follows  ttv.t  the  solutions  of  (1.7)  and  (1.74)  obey  a stability  rela- 
tion with  respect  to  changes  in  the  initial  data  and  inhomogeneous  terms;  if 
ut  + L f(u)  - g(t),  r fc  Cl  , t > 0,  (2. 24a' 

B f(u)  - 0,  x fc  9-0.  , t > 0,  (2.74b' 

and 

u(0)  - u0  , x 4 XL  , (2.24c' 
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,-fc  i 

llu(x.t)  - u(x,t]!\Ll(n)-  II  uo  - u0\lL,ca)  -v-  \ollg(t')  - g(tOllL-CXI)  dt'  • (2.24d) 
Seooud,  using  this  and  properties  (1.10)  and  (1.22)  of  g,  we  see  that  for 
'Y  sufficiently  small  and  mT'  suff iciently  large,  the  solution  U of  (1.24)  may 


be  replaced  with  arbitrary  accuracy  in  L1  C-H-)  by  the  solution  of 

TTt  -V  L f(TT)  g,  ufl  , t > 0,  (2.25a) 

B f(U)  - 0,  x e Cl  , t > 0,  (2.25b) 

m 

u(0)  _ Uc  + Y.  g '(nr)r,  Ifeil  . (2.25c 

n = \ 

Third,  we  use  (1.10),  (1.22),  (1.12b),  the  density  of  d(L)  in  L1  C.D-  ) , 


(2.25c ' 


the  fact  that,  for  u 6 D(L), 
(S(T  ) - I)  u - 


o as  r 


(2.25) 


and  the  L -oontractivene33  of  , to  show  that  for  'f  sufficiently  small 
and  mT"  sufficiently  large,  the  functions  un  of  (1.12)  may  be  replaced  with 
arbitrary  accuracy  in  l_'  Cil)  b y the  functions  un  given  by 


^ o o 
u n u 


(2.27a) 


n-t-t  _ rJ  n 
u - u 4- 


g’((n+l)r)T  n^m 


(2.27C1 


prom  thQ  assumed  convergence  (a 3)  and  the  assumption  that  a steady  3tate 
is  approached,  (a2)  or  (a2*),  we  note  that  the  steady  state  pr  of  the  solution  of 


(2.25)  can  be  approximated  by  the  functions  XT  given  by 


- 


^ nti  ^ n 

u - r , un , 


Y g'(pr)  r,  itil  , 
p=> 


(2.. 23a) 
(2.23d) 


with  arbitrary  accuracy  in  L'  CfL)  if  n T is  sufficiently  large  and  O'  is  suf- 
ficiently small;  similarly  the  steady  state  u of  the  solution  of  (1.7)  can  be 
approximated  with  arbitrary  accuracy  in  UC!!)  for  nY  large  enough  and  O'  small 
enough  by  the  functions  un  generated  in  (2.27). 

rO  n Aj  n 

When  uQ  ^ u and  g*  :j.  0,  (2.23)  shows  that  u and  XT  just  given  satisfy 


min(xjn  ,u+  ) - min.(un  ,u* 


0 as  n — =»  oo  . 


. ,>*  ' ' ' (.H)  * ' 

Finally,  when  we  combine  the  above  observations  with  permrk  2,  we  obtain 


(2.29) 


m 


Theorem  1 easily. 
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3.  Direct  proof  for  the  one-phase  Stefan  problem 


It  would  be  interesting  to  see  to  what  extent  a direct  proof  of  Theorem  1 
can  be  fasnionei,  one  which  has  no  recourse  to  the  algorithm  (1.12)  and  the 
approximate  solution  which  it  generates.  Wo  will  find  it  convenient  to  intro- 


duce, besides  the  dependent  function  u,  the  function  v(x)  defined  by 
v(x)  - f (f(u(x,t))  - h(x))  A-t  , 


(3.1) 


where  h(x)  satisfies  (1.4u,b)  and  (1.9a)  when  appropriate.  From  (1.7),  (1.10) 


and  (Al ) , 


IT- 


(t)dt  +-  UD  - u. 


(3.2) 


From  (1.4c)  or  (1.9b),  v given  by  (3.2),  when  it  exists,  will  be  unique  a.e. 
Theorem  1 says  that  u in  (3.8)  depends  on  uc  and  g'  only  in  the  combination 


U<n, 


when  uD  'z.  u"  and  g*  -5.  0.  Thus,  v also  depends  on  u0  and  g»  only  through 
that  combination. 

Let  us  consider  the  problem  of  detennining  v for  the  one-phase  Stefan 
problem.  This  problem  is  the  special  case  of  (1.7)  with  l -A,  g - 0,  and 


u £ X 


u - X 


(3.3) 


(We  must  have  X 'Z-  0 in  order  to  satisfy  f(0)  m 0.)  yrom  (1.25)  we  have 
u+(x)  m A , u“  (x)  - - 00  , and  thus  the  condition  (2.1a)  is  satisfied  automat- 


ically. 


A typical  set  of  initial  conditions  is  the  following; 

Uo^  A , x €.  JfL0C  n 

ue  = 0,  x 4 /L-.fL0  . 


(3.4a) 

(3.4b) 


Since  we  are  not  looking  for  maximum  generality,  we  cousider  the  case  where  the 


sources  satisfy 


3uppg*(.,t)  c oi(il0)  V t 0.  ( 

The  problem  always  has  a steady  state,  since,  when  g*  ^ 0,  the  set 

n (t)  m mt  ^ x \ u(x,t)^A^ 


(3.4c) 


(3.5a) 
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increases  monotonically  with  t and 


meas(  A (t ) ) ^ -i-(^juodLx  + j'J'g’dxd-t)  V t.  (3.5b) 


The  steady  state  is  of  the  form 

r X 


u 


x e n 

x € intC-O--  -0.) 


where  the  set  -A  satisfies 

A-oafL  cfl 

From  (3.1), 


▼(*)  - f . (u(x,t)  - A ) «lt 


where 

t*(x)  = inf  {%  | u(x,$)-*\]  . 

Clearly 

v(x)  > 0,  x 6 int(  A.  ) , 
v(x)  - 0,  x € fl  - intffl  ), 
and  in  particular 

v (x)  — 0,  x6  3 Cl  - 3 TL 
From  (3.7a)  and  the  Stefan  equation, 


(3.6a) 

(3.6b) 

(3.7a) 

(3.7b) 

(3.9a) 

(3.8b) 

(3.8c) 


I 


3fla)-3A 


- n-V  v dS 


=1  - 

J rzn  - 


- n-V  v dS 


C3rL-3nct>)  n 3A 

+ A meus(  fi  - fl  (t)),  (3.9a) 


where  in  the  first  integral  n points  out  of  A (t),  and  in  the  second  n points 
out  of  A . This  gives  us 

n-7T(i)  ci  0,  x 6 9/1  - 3A  • (3.9b) 

According  to  (3.6),  inside  A v satisfies 


s’(t)  at- 


(3.10) 


A v-  A - uQ  - 

Boundary  conditions  on  v(x)  far  x 6 3 A are  found  from  (1.1b); 

B v - 0,  x « 9-0.  • (3.11) 

If  we  can  show  that  v depends  on  uQ  and  g*  only  through  the  combination 


*T  — u + 

J rs  o 


L 


8*  dt- 


(3.12) 


It  'will  follow  from  (3.2)  tint  u is  also  dependent  only  on  u#  • But  (3.10), 
(3.9 c),  (3.9b),  (3.11),  (3.9u),  and  (3.6b)  define  an  elliptic  free  boundary 


J 
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problem,  and  the  question  becomes  one  of  the  uniqueness  of  its  solution.  This 
elliptic  free  boundary  problem  is  of  a type  which  has  been  studied  in  connection 
with  a diffusion-consumption  problem  (Refs.  1,  8).  with  v expended  to  /X  -D. 
according  to  (3.8b),  the  problem  may  be  written  as  the  determination  of  the 


steady  state  of 


W-t  n A w + TJ0 


w > 0 


w - 0 


w(0)  = w 0. 


(3.13a) 


(3.13b) 


That  is , any  solution  of  the  elliptic  free  boundary  problem  will  also  be  a 
steady  state  solution  of  (3.13)  for  an  appropriate  choice  of  w0  . (Recall  that 
we  have  assumed  supp  TJ  c Cl(  -0-o  )•)  By  regarding  (3.13)  as  the  limit  as 


£ ^ 0 of  a problem  of  the  form 


wefc  m A we  4-  TJ0 


w&  >,  £ 


wfe  x £ 


(3.14) 


one  may  prove  the  monotone  dependence  of  the  steady  stats  on  wD  (Ref.  8).  pn 
particular,  from  (3.3a)  and  (3.3),  the  steady  state  will  satisfy 

V"  x.  v £ v+  (3.15) 

whore  v+  are  the  steady  state  solutions  of  (3.13)  corresponding  respectively 


w0  m 0 


>*=  l 


SI  G(X,X»)  (no  (K*)  - X ) dx1 


A i(x,x*)  m - % (x  - x* ) , x e XI  , x’e-H_  , 

B Gr(X,X»)  m 0,  X 6 dO.  , X*  £ -G.  . 


(Wo  can  show  that 


w“(x,t)  m C f(u(x,t’))  fit' 
Jo 


(3.15a) 


[3.16b) 


[3.16c' 


[3.154: 


( 3. I7) ) 


The  nociotonicity  cited  for  solutions  of  (3.13)  slows  that  the  respective  solu- 
tions w+  corresponding  to  the  initial  data  wQ+  satisfy 

0 £ w“*  - w~  x.  / (3.' 3a : 


(3. '3a: 
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whore  Y satisfies 

rt-  ay 

(3.18b) 

Y (0)  n wo+  , X feil  , 

(3.18c) 

BY  - o,  I (r  2fl  . 

(3.13d) 

By  using  (1.4o)  or  (1.9b),  we  find  that  the  steady  state  of  Y is  0,  and  thus 
from  (3.15) 

▼ " - v - , (3.19) 

that  is,  the  solution  of  the  elliptic  free  boundary  problem  is  unique. 

4.  positive  and  Negative  sources  in  the  One-phase  ste^an  problem 

in  proving  Theorem  1,  we  assumed  that  the  sources  g*  wore  non-negative, 
in  the  case  of  the  one-phase  stefun  problem  where  l -A  , i - 0,  and  f is 
given  by  (3.3),  we  can  show  that  the  theorem  doe3  not  generally  hold  without 
imposing  the  condition  (2.1b).  For  example,  one  may  have  tjo  given  by  (3.12) 
such  that 

U0  \ , * £ Sl0  CL  fL  , (4.1a) 

TT0  - 0,  x£fl-fL0  , (4.1b) 

and  also  g*  such  that 

g*  ^ 0 V x € 0.  and  V t v 0,  J g'dt  < 0,  -H-o  . ( 4. 1c ) 

in  t hi 3 case  the  steady  state  of  the  solution  of  (1.24)  will  have  p = A for 
x A0  • However,  suppose  g»(x,t)n  0 for  t <:  Tc  . If  TQ  is  large  enough, 
the  solution  of  (1.7)  will  approach,  as  t — > to  , a function  which  is c A+  £ 
everywhere,  where  t > 0 and  € — =»  0 as  T0 — * oo  . Then  the  addition  of  the 
sources  gf  for  T0  < t <■  oo  will  result  finally  in  a steady  state  u wnich 
satisfies  u-uA  for  x e .Q.  j . Clearly  tnis  is  not  the  same  as  TT. 

The  proof  given  in  the  last  section  of  the  uniqueness  of  v satisfying 
(3.11),  (3.8c),  (3.9b),  (3.11),  (3.8a),  and  (3.6b)  remains  valid  even  in  the 
cane  of  negative  sources.  However,  when  the  requirement  g*  ^ 0 is  dropped,  the 
monotonia ity  stressed  by  (3.5a)  no  longer  holds,  and  the  steady  3tate  u need 
not  be  given  by  (3.6a'.  accordingly,  the  equation  (3.2)  satisfied  by  v will  not 
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load  to  (3.10)  In  general,  but  only  when  u-  A where  v > 0.  since  v 7 0 at  a 
point  x means  that  at  some  time  t,  u(x,t)  , we  can  insure  that  u vd.ll 
satisfy  this  requirement  if  for  all  x t IT  and  t > 0 we  have 

> 0.  (4 .3) 


More  generally,  let  u and  u be  solutions  of  the  one-phase  Stefan  problem 
-with  initial  data  uG  , u0  and  sources  g*,  g* , respectively,  such  that 


= u0+  i g.at  =.  t 

b J n 


g' 


(4.3) 


and  let 

n°° 

< 

il 

o' 

f(u(x,t))d-t  , 

▼=  J0  f(u(x,t))dt:  . 

(4.4' 

Define 

<5-=  ^xfe-Q. 

r°°  1 

l v(x)  > 0 and 

j-t  g*  (x,t  * )dt  0 for  some  t > 0 J , 

(4.3a' 

cr  5T  \_xt  -0- 

| v(x)  > 0 and 

g'  (x,t*  )d‘t/  0 for  some  t ol  . 

(4.5b) 

If 

u(x)  - u(x) 

for  all  1 fe  c U a , 

(4.5a) 

then 

u'x)  - 

u(x)  for  all  x fe -Q.  . 

(4.5o) 

This  result  is  derived  easily  from  the  following  observations. 

First , 

if  x (=.  fi  - o’ 

f 00 

, either  v(x)  - 

0 or  v(x)  -?■  0 and  Jfc  g*(x,t*)  dt'  ~2- 

C for 

all  t ^ 0.  in 

the  latter  case 

, u(x,t)  >A  for  seme  t > 0,  and  the 

c or.iltlon 

ou  gt  implies  that  u(x)-  A . 

Seco.i,  v satisfies  the  following  elliptic  free  bound  . ry  . ro'.  'er.; 

v(x)  > 0 for  x t int'fl  1,  (4.,'a' 

3 [*|u0  + 10  6*  it?- A]  , 

v(x)  ■=  0,  x t 3 A - 9 i).  , 4."e ' 

n • V v ( x ) n 0,  xfc  9/1  - ?fl  , (4  -M' 

E v - 0,  x 6 3-0.  , (4.  ’’a' 

A v - u - T70  1 xfe  '.  (4.nf' 

Since  wo  only  compare  solutions  of  f 4 . "*  * f v ir.  7.  t es  » re- 

prescribed  value  at  3 jne  points  o'*  A and  -.s  \ ver  t..«  re-'.  1 , • e 


uniqueness  of  v follows  as  before.  By  '3.1'.  t -»  u.i^.e  .vs 


s oat  >- 
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lished  under  these  conditions.  Note  that  we  no  longer  have  the  restrictions 
(3.4)  on  uQ  and  g». 

It  is  also  possible  to  verify  these  facts  by  direct  consideration  of  the 
approximating  algorithm  (1.12)  and  invocation  of  the  convergence  of  the  approx- 
imate solution  to  the  exact  solution. 


5.  Application  to  a problan  of  Anodic  Smoothing 


Consider  the  Stefan  problem 
u-t  = A f(u) 

whers  f (u)  is  given  by  (3.3)  and  with  initial  data 
uD  n A -t-  dicx,  £"0  , x d n 

uQ  - 0,  x e fL-fi0 


(5.1) 


(5.2a) 

(5.2b) 


where  £ 0 and  do.  is  a small  positive  number.  Let  ns  impose  the  boundary 


condition 


f(u)  - 0,  x t 3-0  . 


For  this  problem  we  use  the  notation 


= £ f(u(x,t))<tt  . 


(3.10)  becomes 


A (dv)m  — do«.  %0  , x 0.o  , 

and  to  lowest  order  in  d<*  , from  (3. Bo)  and  (3.9b), 


dv  no.  0,  ife  3-Qq  — 3-0. 


In  addition, 


dv  -m  0 , x 3 A.  . 


Thun , 


n-V  (dv)=dk  J n-V  0o  (x,x»)  f0(x*)dx'  . 


wl-ors  is  given  by 


0o  (x,x*)  - " S(*  - *’)  . x 6 -O-  0 ,x»fe/l0  , 
0 (x,x* ) r 0,  xe  3-0.o  H 90  , x*t  n. 


(5.2c) 


(5.3) 


(5.4a) 


(5.4b) 


[5.4c' 


(5.5- 


( 5 . > 


, x’fe  fl0  , (5.6b) 


0o  (x,c»)  - 0,  x 6 - 3/1  , x’feHo  • 

denote  the  steady  state  domain  A-  by 


(5. del 
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TT-  - IX  . (5.7) 

From  (3.9b)  and  (3.10),  (5.5)  Implies  that  Udot  is  obtained  from  IT0  5/ 
traveling  outward  normal  to  IT0  at  each  point  of  9 IT 0-9H.  a distance 
-A_  d*.  ( n-VG0  (x,x*)  (x')dx' 

* JH0 

We  may  now  envisage  a collection  of  such  problems,  in  which  we  generate 
domains  jQ  «.  IT  , aL  ■>  0,  where  H <*_-*. «-  is  found  from  .fi. 

by  traveling  outward  along  tne  normal  to  9 -0.^  at  each  point  of  911^-9-0. 
a distance 

_ -J-  ^ a-7  = V(oO<Aot  , (5.3a) 

where 

A v^  - “I’o^  . x 6 IT*  . (5.8b) 

v^  = 0,  it  9-11*  , (5.8c) 

and 

f ^ 0,  x fe  -0.  * , (5.9a) 

^:0,  it  il-H  * . (5.9b) 

This  collection  of  problems  may  be  tnought  of  as  defining  a ’’time-dependent” 
free  boundary  problem,  where  ot  represents  the  "time”  variable. 

The  re3ilts  of  sections  3 and  4 show  that  the  region  IT  ^ can  be  obtained 
directly  as  the  support  of  the  steady  state  u of  the  solution  of  (5.1)  with 
initial  data 

a. 

uQ  - A -*-]q  da.1  , xt  H o , (5.  iDa) 

u o - r*  £.<  <*<*.'  rxeiT-XI0  . (5.19b) 

The  "time-dependent"  free  boundary  problem  Just  described  is  very  similar 
to  a problem  which  occurs  in  the  theory  of  anodic  smooth j ig  (Refi.  3,  5,  6,  7). 


To  bring  problem  (5.8)  inti  the  form  of  the  anodic  smoothing  problem,  we  find 


a function  ©* 

witn  the  following  properties; 

(x)  n 1,  x fe  9 n A 9-ft*. 

(5.11a) 

(X)  X 0,  x 6 IT-  IT ^ 

(5.11b) 

bounded  and  non-negative,  x fe  IT  . 

(5.11c) 

Next  let 
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Wlt.i  given  by  (5.*d),  the  function 

%.  - V + 6<* 


has  the  following  properties; 


(5.12) 


(5.13) 


z 0,  x fcJl*  , (5.14a) 

- 1,  * 6 9/1  rt  ^-0-^  , (5.14b) 

0,  i€  9X1*  "3X1  , (5.14o) 

and  XI A is  generated  by  traveling  outward  along  the  normal  to  9X1^/  , for 

0 £ <Z  <*  , at  eaoh  point  of  3X1^/  -9X1  , at  the  rate  (with  respect  toot'  ) 


v - "T  n'v  ^ * 


(5. 14d ) 


in  the  anodic  smoothing  problem  (Ref.  3),  XI  is  the  domain  exterior  to 
a cathode  c,  3X10  fi  3X1  ■=.  9_fl  j*  H 3/1  ~ £c,  for  o^-^,  0,  and  A-  1.  Thus 
we  3ee  from  (5.11)  th-it  ©0  may  also  be  used  as  ©<*  , 0,  and  from  (5.12) 

I-  may  be  used  as  for  <*  > 0.  in  that  ease  the  initial  data  for  an 


equivalent  steady  st^te  stefuri  problem  are 

Uq  A + OC  I X 4 XI  Q f 

U - — 0,  X €:  XI  — XX  q . 


(5.15a) 

(5.15b) 


It  al30  follows  that  any  effective  methods  for  solving  the  steady  state 


Stef  an  problem  (5.1)  and  (5.15)  (Ref.  2)  can  ba  used  to  solve  the  anodic 
smoothing  problem,  in  addition,  regularity  results  for  the  free  surface  in  the 
steady  state  Stefan  problem  carry  error  immediately  to  the  fre?  surface  for  the 
anodic  smoothing  problem,  independently,  Charles  V.  Elliott  has  made  similar 
observations  about  the  equivalence  of  the  anodic  smoothing  problem  and  a steady 
state  stefan  problem  (Ref.  4). 
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